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Abstract 

Conductance through a system consisting of a wire with side-attached quantum dots 
is calculated. Such geometry of the device allows to study the coexistence of quantum 
interference, electron correlations and their influence on conductance. We underline 
the differences between "classical" Fano resonance in which the resonant channel 
is of single-particle nature and "many-body" Fano resonance with the resonant 
channel formed by Kondo effect. The influence of electron-electron interactions on 
the Fano resonance shape is also analyzed. 
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1 Introduction 



A model consisting of a quantum dot (QD) attached to a metallic quantum 
wire is conceptually analogous to the Fano model [1] consisting of a contin- 
uous spectrum and a discrete level. Such geometry is also in close relation 
to the behavior of magnetic impurities embedded in the host metal, where 
the conduction electrons are scattered resonantly on impurities while propa- 
gating through the metal. This arrangement can be opposed to the standard 
geometry, where a quantum dot is connected in series with leads and reso- 
nantly enhances the conductance. A model of a quantum dot side-coupled to 
the quantum wire is also applicable for the simplest description of a mag- 
netic atom deposited on the metallic surface. In such a "mapping" the QD 
in Kondo regime mimics the adatom with nonzero magnetic moment and the 
wire (when assumed to be multichannel) serves as an analog of the metal- 
lic surface host. Thus, the model considered by us can provide information 
applicable to various nanostructures. 
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Fano effect, described in 1961, has been observed in many physical systems: in 
nuclear scattering [2], scattering from donor impurities in an electron waveg- 
uide [3], optical absorbtion [4], electronic transport in heteroj unctions [5]. 
Recently it is seen in new light due to the rapid progress in nanotechnology. 
Apart from 'classical' Fano effect, where the resonant channel is provided by 
discrete single-particle level, a new type emerged, where the resonant channel 
is formed by many-body effects. This additional resonant channel is formed 
by the electron-electron interactions resulting in the sharp Abrikosov-Suhl 
(Kondo) resonance in the vicinity of Fermi surface. This kind of Fano resonance 
has been observed recently in scanning tunnelling microscope experiments for 
single atoms placed on metallic surface [6,7] and also in lateral quantum dots 
[8]. There were theoretical attempts to describe these phenomena in [9] and 
in [10], respectively. 

A single quantum dot modelled by Anderson Hamiltonian and side-attached to 
a perfect wire was investigated in resonant regime by slave boson mean field 
approach [11], exact diagonalization [12], interpolative perturbative scheme 
(IPS) [13] and X- boson treatment [14]. For this arrangement we concentrate 
on the crossover between "classical" and "many-body" Fano resonances. 

A system of two quantum dots placed in arms of Aharonov-Bohm ring has been 
investigated recently for non-interacting electrons [15]. In our model we study 
the influence of the strong electron correlations within the dots on transport. 
Repulsive Coulomb interactions in the nanoscale devices are defined by the 
charging energy ~ e 2 /2C, where C is the capacitance of the device and e- 
electronic charge. They cause the most celebrated phenomena as Coulomb 
blockade and resonant (Kondo) electron transport [8]. 



2 Model and calculations 

Hamiltonian of the nanodevice consisting of metallic wire and quantum dots 
attached parallel to it, is taken in the form: 

H = E € k4,a C k,a + EE e 7<<A- 

fc,cr 7=1,2 cr 

+ E U ~i n ^ n ~tl + E t 7 [d+ a c ka + h.c] (1) 

7=1,2 7=1,2 

cr,k 

Creation operator c^ a {d^ a ) describes the conduction electron in the wire with 
momentum k, energy and spin a (7-th QD localized electron with spin a). 
n 7)(J is the number of electrons with the spin a at the localized state 7 with 
the Coulomb repulsion U 1 . The last term describes hopping between QDs and 
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wire. We assume spin-only degeneracy of the discrete QD's levels and the wire 
to be single-channel. The subscript 7 represents the orbital quantum number, 
which follows from quantization due to spatial confinement of the 2D electron 
gas within the quantum dot (see for example [16]). We underline that 7 is not 
conserved in the process of hopping of electron between QD and causes an 
indirect interaction between the dots. 

To investigate transport through nanodevice, the wire has been connected 
to external electrodes: H con = ^ a =L,R t a [c^ aa Ck' a + h.c.]; c£ aa is the creation 
operator of the electron with the spin a and energy e^ in the electrode a. 
Spectral densities of electronic states in the wire p w o i(T (u) and in the electrodes 
p e i^{u) have been assumed to have Lorentzian shape with a halfwidth much 
larger than Kondo temperature of each QD. 

The numerical parameters of the equal Coulomb repulsion U 7 — 1 meV (for 
7 =1, 2) and tunneling strength T 1)(T (uj = e F ) = 2itt 2 1 p w , a (uj = e F ) = T imax ^ = 
0.28 meV have been chosen to meet the experimental QD's data (e F = 
defines the Fermi level in the wire). A scheme of the investigated nanodevice 
is presented in the a-inset of Fig. (2). 

Following Meir and Wingreen [23], conductance can be related to the spectral 
density of the central region (in the present case - metallic quantum wire with 
quantum dots attached to it) and its coupling r L ( R ) ja (u) = 27it 2 L ^p e i^(u) to 
left (L) and right (R) electrode. For the symmetric coupling and in the zero 
limit of drain-source voltage conductance can be written in the form: 

S = ^ E / IV(e) (-^f) Pw Ae)de, (2) 

a —00 ^ ' 

where /(e) is the Fermi distribution function and r CT (e) = T Lj0 .(e)T R (T (e)/(T Lja .(e) + 
r_R,(r(e)). Spectral density of the investigated nanodevice is calculated from 
the appropriate retarded Green function: p w ^{uj) = — (l/7r)ImG CT (u; + i5). At 
T = 0, when —^7^ — > 5(e) the formula for conductance takes a simple form 

S = ^£,r CT (oK, CT (o). 

To calculate conductance through the considered nanodevice the Green func- 
tion of the electronic wave propagating through the wire should be calculated 
in presence of QDs attached to the wire. We start from the Dyson equation 
(energy and spin dependence of the Green functions and selfenergies is under- 
stood): 

G = G + G SG, (3) 
which can be written by iterating in the form of T-matrix: 
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G = Go + GoTGo, T = E/(l - £G„). 



(4) 



Go is a free conduction electron Green function. The form of selfenergy (and 
T-matrix) depends on the geometry of nanodevice and approximation made. 

Dyson equation describing scattering of the electronic wave on both quantum 
dots: 



G — Gq + Gq^iGq + GoS 2 Go + Gq^iGq^Gq 

+GqTj2GqT,iGq + ... (5) 

can be cast into the form Eq.( 3), when £ = Si + E 2 . The selfenergy E 7 , (7 = 
1, 2) arises from scattering of the conduction electron propagating through the 
wire on the 7-th quantum dot. Thus, the total selfenergy £ contains infinite 
number of scattering events in quantum dot 7 alone and multiple scatterings 
involving both quantum dots. Utilizing the second part of Eq.(4) we can ex- 
press £1 and £ 2 by T x and T 2 and the selfenergy due to both quantum dots 
is: 

S = ? S7= (l+T 1 Go) + (l + T 2 G )- (6) 



For non-interacting electrons, the expression for T-matrix due to one QD 
attached to the wire obtained by equation of motion (EOM) for conduction 
electron Green function is of the form T 7iCr (o;) = t^GQD-y^i^), where the Green 
function of the dot Gq^ 7 (w) represents a simple single-particle level. In the 
interacting case, it is replaced by dressed Green function of the localized QD 
level, which has to be calculated in such approximation that properly describes 
Kondo effect. In the present work the Interpolative Perturbative Scheme (IPS) 
is utilized [17,18,19]. It is based on the selfconsistent second order perturba- 
tion (SOPT) in Coulomb repulsion U [20,21]. However, within the IPS the 
calculated selfenergy due to strong electron correlations E^ 5 is taken in the 
form, that has correct limits: when U — > 0, £ /P5 — > Y> SOPT and when r — > 0, 
T, IPS — > Y] atom , This allows to calculate the conductance both in resonant 
regime and in Coulomb blockade regime. Moreover, by replacing Hartree-Fock 
solution for the impurity level by an effective field determined selfconsistently 
from the condition on particle number, the fulfillment of the Friedel-Langreth 
sum rule [22] is also obtained within this method. It is applicable both for 
T = and finite temperatures. 

If the expressions for T x and T 2 -matrices are known, then the total selfenergy 
Eq.(6) is calculated and finally total T-matrix by second part of Eq.(4): 

rp = T 1 (1 + T 2 G )+T 2 (1 + T 1 G ) 

' (1 + T 1 G )(1+T 2 G )-G [T 1 (1 + T 2 G )+T 2 (1 + T 1 G )Y U 
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The extension of the above formalism on 7 = iV QDs connected to the wire 
in a star-like fashion is straightforward. 



3 One quantum dot attached to the wire 

Utilizing Eq.(4) spectral density for the propagating electron can can be writ- 
ten in the following form: 

PwA u ) = PwoA^) x 
{1 + ImG 0><T (u)[ImT ff (u)(^ - 1) - 2q a ReT a (u;)}}, (8) 

where we have denoted: 

ReG 0:a (u + i5) _ 

qa= ImG 0ttT (u + i8)~ r ff ' U 

which can be identified as Fano parameter. The real and imaginary parts 
of the free retarded conduction electron Green function are: ReGo, a (ou) = 
P J dep w0tCT (e) I (uj — e) and ImG 0ja (ui) = —np^^^u), and parameters A CT = 
t^-ReG ,<T and the r 7j(7 = id"^p W Q^{uj). In the non-interacting case when the QD 
is represented by a single-particle level, Gqd^o- — [oj — e 7)(7 — A CT + iY^^Y 1 , 
and making substitution e = {u — e 7j(T — A cr )/r 7j(T we get the Fano well-known 
formula from Eq.( 8): p w>er (u) = p w0>er (u})[(e + g CT )V( e2 + 

For the interacting quantum dot, the Green function in expression for T-matrix 
should contain appropriate information on correlations. Within IPS it has the 
form: 

= [co - e 1>a A CT n^U - S(f + ir^]' 1 . (10) 

Thus, the bare quantum dot level is additionally renormalized and gets the 
width due to selfenergy S /P5 by electron-electron interactions and the Kondo 
peak is generated by specific cu-dependence of S /P5 . 

Numerical calculations of the selfenergy within IPS have been performed under 
assumption of smooth and weakly u;-dependent density of states in the wire 
near the Fermi level. This is a requirement for stability of numerical procedures 
when calculating selfenergy S /PS '. Then, the integral defining ReG 0i(T (uj) = 

D 

P I dep w0:a (e) I \uo - e) = p w0i Jn\(D + uo)/(D - uo)\ -> for halfwidth of 

-D 

spectral density of the wire D 3> uo implying q a = (see Eq.(9)). 
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Fig. 1. Conductance through the wire with one QD side attached to it vs. QD spa- 
tial level position. Two regimes are presented: resonant regime (RR) and Coulomb 
blockade (CB) regime. Calculations are made for Coulomb repulsion U =lmeV. 
Open circles - U =0 case. 

Utilizing Eqs. (8) and (2), at T = (and for q a = 0) the formula for conduc- 
tance can be written: 



27re 2 
h 



i - 



r 2 

lmax,a 



V c cr ) ~ x lmax,<r 



where for non-interacting case e CT = ei j(T is the bare QD level, and in presence 
of interactions e CT = e e ffi,a = e i,o- + n\- a U + -ReS{ PS (0) is the renormalized 
level. In this particular case, when IrriE IP (u = 0) =0 as for Fermi liquid, 
the conductance for non-interacting and interacting cases can be expressed by 
the same relation. It is a consequence of the fact, that Hamiltonian Eq.(l) in 
this case can be formally written as for non-interacting resonant level model, 
with bare QD levels e 7)(7 replaced by renormalized ones e e // 7)(7 . Conductance 
vs. QD level position is plotted in Fig.(l). Two types of curves are shown: in 
resonant regime (RR), with one broad antiresonance, and in Coulomb blockade 
(CB) regime with two antiresonances. In the resonant regime, for T = 0, the 
total extinction of conductance is observed for e 7 < and e 7 + U > 0, i.e. 
when the QD level is singly occupied. In this case the destructive interference 
has its full strength due to Kondo resonance. For e 7 > (empty level) or 
e 7 + U < (doubly occupied level) the dot weakly disturbs transport through 
the wire and conductance approaches unitary limit. For this configuration 
Kondo temperature T K ~ 700 mK for e 7iCT = —U/2, as obtained from the 
full width of the QD density of states at half maximum. When temperature 
increases, the Kondo effect is gradually diminished and the scattering of the 
electronic waves by many-body QD's Kondo level becomes less effective. By 
a decrease of QD coupling to the wire the system is driven into the Coulomb 
blockade regime. By this operation the Fano resonance caused by effective 
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Kondo peak in the vicinity of Fermi surface is depressed. Instead, two Fano 
resonances appear due to scattering on two single particle levels with the 
separation of the order of U. Although the Kondo temperature of the QD 
has been decreased by lowering hopping t±, the electron correlations still exist 
in the system. When the temperature is increased the CB peaks get sharper 
due to further diminishing of strong electron correlations. For comparison, the 
conductance for U — has been calculated. In this case the total suppression 
of conductance takes place when the bare dot level crosses Fermi level. 

How Fano resonances appear conductance in Fig. (1) can be better understood 
in a simplified picture. At temperatures much lower than Kondo temperature, 
T <C T K , the Green function of interacting quantum dot can be approximated 
as: 



G'M > + iS) = \ — 

+ + ^ (12) 

uj - e ha -A a -U + iT ltff uj-e K + iT K 

The first and second terms simulate charge peaks in the spectral density and 
the third term represents many-body Kondo peak of the width proportional 
to the Kondo temperature Tk and position ex in the vicinity of Fermi surface. 
Zi, Zjj and Zk are appropriate strengths of the poles: Z\ + Zjj + Zk < 1. 
In the Kondo limit the Kondo peak can be regarded as a potential scatterer 
[24] and Z K ~ irT K /T [25]. Utilizing Eqs.(12) and (8) the spectral density of 
conduction electron scattered by QD is written: 



/ x / x (, . 7 Ql -2q aXl -l q 2 a - 2q a x K - 1 \ 
PwAu) = PwoAu) I 1 + Zx -z—^ + 7T x 2 + i )' ( 13 ) 

where following substitutions were made: x\ = w ~ e ^~ a and xk = ■ The 
second and third terms produce the shape of Fano resonance. When temper- 
ature is sufficiently low and coupling to the wire is large, electron interactions 
prevail in the dot. Then the third term in Eq.(13) dominates and the Fano 
resonance in the conductance vs. position of the QD spatial level dependence 
is due to many-body effects. The maximum of Kondo resonance (which is 
transformed into Fano resonance in conductance for present geometry) takes 
place for e\ = —U /2, i.e. for symmetric Anderson model. The influence of the 
second term describing "classical" Fano effect is practically obscured by the 
Kondo effect. It becomes to be visible in the Coulomb Blockade regime, when 
Kondo temperature of the system is decreased by lowering the hybridization 
strength T 7 of QD with the wire. Thus, the Kondo temperature is the energy 
scale which governs the transition between "classical" and "many-body" Fano 
resonances. 
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Phase shift acquired due to scattering on the localized level can be derived 
when the form of T-matrix is known [24]: i](uj) = argT(u + i8). For non- 
interacting electrons, at the Fermi level ou — 0, it has the form r)(ui = 0) = 
arctan(j^-), where e^ a is the bare QD level. Thus, the phase is changed by 
tt/2 when QD level e lj(7 crosses the Fermi level. For interacting case rf(u = 0) = 
arctan(^), because the resonant level is of many-body nature. The position 
of the Kondo peak is 6k — 0, which gives r](oo = 0) = 7r/2. 



4 Two quantum dots side attached to the wire 



The formula for the Green function of conduction electron scattered by two 
quantum dots was obtained from Eq.(4), for the selfenergy given by Eq.(6) 
and total T-matrix given by Eq.(7). The spectral density has the form of: 



PwA u ) = PwoA^) x 

[(U - € ha ) (u - e 2 , a )} 2 + [ r i,a(^ - €2,a ) + ^2,a{ U ~ 

and conductance at T=0 can be expressed as follows: 



27re 2 
h 



E r -( )p-o, CT (o) 



1 - 



(14) 



[^l,a^2max,a ^2,cr^- \max ,<t] 
(ei,o-e2,o-) 2 + [^l,a^2max,a + ^2,a^lmax,aY 



.(15) 



Thus, Q = when e± = or e 2 = is tuned to the Fermi surface in the wire. 
On the other hand, for e 7 = (— r 7 /Ty)ey conductance reaches unitary limit at 
T = 0. For such condition the total phase shift r)(u) = —arctan(J2 1= i,2 
is equal to zero and the electronic wave propagates through the wire without 
scattering. Thus, by setting a proper configuration of the QD's levels by gate 
voltages, the conductance through device can be set to unity or to zero. 

When the level of one QD, say e 2 , is kept fixed, while e\ is shifted by the 
gate voltage, then the Fano-like resonances appear in the conductance. The 
QD level e 2 acts as additional background channel and the shape of Fano 
resonance depends on its position with respect to the Fermi level. Indeed, 
equation (15) can be approximated by the expression: 



2vre 2 

6~ T ^r„(o)^(o) 

where following substitutions have been made: e a = e 1)Cr /T ma:El>cr and Fano 
parameter q a = e 2) „lV max2)a . 



r 



2 

max* 



(e CT + qaf 



~2,a 



4- V 2 

' ^ max2,a 



(16) 
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The non-conservation of the orbital quantum number while hopping of elec- 
trons between quantum dots and wire is crucial for obtaining Fano-like shaped 
conductance. In the opposite case no information could flow from one quantum 
dot to the other and the dots would act as separate entities. 

To study this problem in detail we considered the Hamiltonian, Eq.(l), without 
Coulomb interactions which is quadratic and can be solved exactly. EOM 
method gives the expression for Fourier transformed retarded Green function 
G 7 , a (t-t>) = ({d 7 , a (t);d+ a (t'))) of 7-QD level: 

= — , T 1 , r 7 , CT r y , CT - (17) 

to e T a -+- 11 7j(7 -+- iV 



Thus, apart from the broadening caused by the hopping of electrons between 
7— QD and wire, the 7-level is additionally shifted and broadened due to 
indirect interaction with 7'-QD level. 

Fourier transformed Green function G^i^t — t') = ((d 7>a (t)] dy a (t'))) describ- 
ing virtual hopping between 7 and 7' levels has the form: 

Gyy C r(^- 7 ) = 7 77 x ft; 7 — : x 7; 7 77? (18) 

{u - e 7j(T )(w - ey ;CT ) + «[r 7i(T (cj - e 7V ) + T^ >(7 (uj - e T>CT )J 



where r 77 / )(T = irt*ty p w o, a - Utilizing Eqs. (17-18) and assuming t 1 (7 = 1,2) 
real, the expression for the total T-matrix due to both quantum dots can be 
rewritten as: 



T a — t\GQDl,a + t\GQD2,a + tit 2 Gi2 t a + tl^\G 2 \^ 

= T 1)CT + T 2i(7 + 2T 12 , CT , (19) 

where T 77 / = t^tyG^y^ and T 77 / j(T = T 7 / 7 <T . The spectral density, Eq.(14), can 
now be written in a more transparent form: 



pwA u ) = PwoA^) I 1 + F ha ImG QD1 (u}) 

+ T 2 ^ImG QD2 {u) + 2T 12 ImG 12 ] . (20) 

Eq.(19) consists of terms describing scattering on each 7 quantum dot sepa- 
rately (but with renormalization due to indirect interaction with 7' quantum 
dot), and the term of multiple scattering. When the 7- quantum dot is gradu- 
ally decoupled from the system (by appropriate decrease of coupling t 7 , or by 
shifting e 7 far from the Fermi level) both the 7-term and last multiple scatter- 
ing term vanish and Eq.(20) coincides with Eq.(8) (with q a = 0) for a single 
quantum dot. When both 7 = 1,2 levels are close to the Fermi surface, the 
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Fig. 2. Conductance through a wire with two QDs side- attached to it vs. spatial 
level position e\ in QD\ for various values of level €2 in QD2. Inset a - scheme of 
the device, inset b - dependence of renormalized level e e //2 level vs. bare €2 level 
position. 

electronic wave traveling through the wire scatters resonantly in both quantum 
dots. Classically, this process can be seen as a multiple bouncing of the ball 
between two walls. This process is summed up to infinity by the Dyson equa- 
tion Eq.(5). The maximum of the "bouncing effect" on the Fermi surface (for 
lo — 0) due to the last term in Eq.(20) takes place for e 7 = rp + p 7 2 y /2 

i.e. when the levels lie below and above Fermi surface. 

When strong electron correlations within the dots are included, the explicit ex- 
pression for conduction electron Green function Eq.(4) has complicated form, 
but it simplifies considerably at T = (see the comment below Eq.(ll)). 
The conductance can then be written in exactly the same form as for non- 
interacting case (Eq.(14, 15 and 16)), where single-particle e 7iCT levels have to 
be replaced by renormalized ones: e e // 7)t7 = e 7jCr + n 7) _ CT C/ + i?e£ 7 P5 (c<j = 0) 
for 7 = 1,2. The discussion of Eqs.(19 and 20) also holds with the above 
replacement. In this special case when IrriL IPS = and the bare QD levels 
are renormalized by static field H e jf = n^- a U + ReYJ PS {uo = 0), all the dy- 
namical many-body effects are erased. In general, apart from renormalized by 
interactions single-particle levels, Kondo resonances are present in the vicinity 
of Fermi surface which have dominant influence on conductance. 

Conductance dependence on the e\ QD level position for various values of level 
e 2 is plotted in Fig. (2) at T = and in presence of correlations. The change 
of the Fano resonance shape can be observed as e 2 takes different values. A 
mirror reflection of the curve takes place when e e //2,(r changes sign changing 
q a = e e ff2,a/T2,a parameter. The e e //2 dependence on bare e 2 is shown in the 
fo-inset. 
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[meV] 



Fig. 3. Conductance through a wire with two QDs side-attached to it vs. spatial level 
position ei in QD\ for a fixed value of level e2=-0.4 meV in QD2 and U =1 meV 
within the dots. Apart from T=0 curve, two curves are plotted for temperatures 
much higher than the Kondo temperature. Inset: conductance curves for interacting 
and non-interacting cases calculated for the same Fano parameter q. 



To compare the conductance curves calculated with and without interactions 
we show two of them in the inset of Fig. (3). The curve for U — has been 
plotted for e 2 equal to corresponding e e //2 found selfconsistently in interacting 
case and the same F2 jC r to have the same Fano parameter in both cases. Two 
substantial differences are seen between the curves: i) the curves are shifted; 
this is caused by renormalization of bare e\ level by interactions, ii) an addi- 
tional plateau emerges for U 7^ which is caused by Kondo resonant scattering 
when the e\ level is occupied by one electron. 

The q = case deserves a comment. In this case e e //2 coincides with Fermi 
level, which corresponds to the bare e 2iC r = —U/2 i.e. maximum Kondo reso- 
nance in QD 2 . For such level arrangement conductance Q = for any value 
of shifted e\ because of full destructive interference caused by QD 2 . This is 
a different situation comparing to classical Fano expression which for q = 
describes a symmetric dip. This shows a resonant nature of the background 
channel provided by the second quantum dot. 

In Fig. (3) we show conductance curves for fixed QD 2 level position e 2 =-0.4 
meV at T = and at temperatures much higher than Kondo temperature. At 
elevated temperatures the plateau in conductance due to Kondo effect on QD 1 
disappears, but the shape of the Fano resonance practically does not change. It 
leads to the conclusion that the shape of the Fano resonance is determined by 
the renormalized single particle level e e //2 and not by the many-body Kondo 
peak caused by resonant scattering on QD 2 . 
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